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1 Introduction 



In non-relativistic mechanics the system or its state can be characterized by 
a number of quantities with dimension of time. Some, such as a resonance 
hfetime, are state independent, and derived from the Hamiltonian; other 
characteristic times depend on the initial state and its dynamical evolution. 
Examples are the dwell, traversal or arrival times, and generally the times 
required to find a dynamical variable with a given value . The quantization 
of the classical characteristic times is problematic and has attracted consid- 
erable attention in recent years |Q. "Traversal" and "arrival" times have in 
particular been studied with many different formal treatments. We shall be 
mainly concerned with "arrival times" in this paper although the difficulties, 
techniques and ideas can be easily translated to other characteristic times. 

1.1 Arrival times in classical mechanics 

Assume a structureless particle that moves in one dimension with position q 
and momentum p. In simple cases the trajectory may cross a given spatial 
point X only once 0, but if a potential barrier reflects the trajectory, or if 
time dependent or noisy interactions affect the motion, the trajectory may 
cross X several times. The first crossing determines the first passage time, 
and in general, the n-th crossing gives the the n-th passage time. For sta- 
tistical ensembles of noninteracting particles there is a distribution of times 
associated with the n-th passage. Provided some restrictions are imposed on 
the initial ensemble, the average value and higher statistical moments can be 
defined from it. However, the distributions may exist even if the moments 
do not. 

Let us examine first the free motion case. An ensemble of freely mov- 
ing noninteracting particles will be described by a phase space distribution 
F{q,p,t), normalized to one, that satisfies F{q,p < 0) = 0, i.e., all par- 
ticles move rightwards. For free motion the trajectories cross the point X 
only once, and J{X, t), the probability flux or current density at X, provides 
the distribution of first passage or arrival times; J{X,t)dt is the fraction of 
particles that cross X between t and t + dt. Using 



(all integrals go from — oo to oo unless indicated otherwise. In this particular 
one the lower limit can of course be set to zero), and the trajectory equation, 
?(^) = Qo + pt/^j the average arrival time is given for free motion by 




(1) 



I J{X,t)tdt = llFiqo,p,0) 



(X - qo)m 



dqodp . 



(2) 



P 
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The integral exists if the singularity at p = is cancelled, at last partially, 
by F. An asymptotic behaviour F ~ p*^, e > 0, when p — >■ is sufficient. 

When the ensemble contains particles with momenta of arbitrary sign, or 
for more complex dynamics with multiple crossings, J is not the ffist passage 
time distribution any more. Suppose however that all particles are initially 
on the left side and that any particle that crosses X is eliminated at the 
crossing instant. When this "absorbing boundary" is imposed it is still true 
that the flux Jab{X,t) = —dN/dt provides the (unnormalized) ffist passage 
time distribution. Here, = N{t) is the time dependent (diminishing) 
"norm" N = JJ F{q,p,t)dqdp, and 

JabiX,t)=VimJ{X-e,t), e>0. (3) 

(Equivalently, the particles that have crossed X at least once may be la- 
beled, without affecting their dynamics, so that the ffist passage time distri- 
bution becomes proportional to the flux of the complementary subensemble 
of unlabeled particles.) If not all particles are eventually absorbed, the dis- 
tribution is normalized by dividing Jab{X,t) by the total norm absorbed, 
l-Ar(oo) =/dtJ„,(X,t). 

A standard way to study ffist passage time distributions in classical 
stochastic diffusive processes is the use of "renewal equations" P]. Assume 
that all members of the ensemble have initial position qo = xq sd. t = 0. 
There are two ways to arrive at point x < X: Directly, via trajectories that 
remain in the left subspace (with respect to X), or crossing X one or more 
times. Let Po{x,xo,t) be the probability density for the "direct" event. It 
can be obtained by solving the diffusion equation with absorbing boundary 
conditions at X. The probability of the "indirect" event can be decomposed 
according to the ffist passage time at X so that the probability density for 
being at x at time t irrespective of the path may be written as 

P{x,xo,t)=Po{x,xo,t)+ f fx,xo{t')P{x,X,t-t')dt' . (4) 

Jo 

(A related equation was proposed by Schrodinger [^,0 in terms of cumulative 
probabilities.) This is the renewal equation, where fx,xo (f) is the ffist passage 
time distribution, and P{x, X,t — t') is the probability density for being at 
X at time t conditioned to having been at X at time t' |]^. fx,xo{t') is the 
probability flux at X with an absorbing boundary, see (||) , but it can also be 
obtained solving the integral equation by Laplace transform. 

1.2 Quantum arrival times? 

Due to the basic role of the trajectory concept to define "arrival" or "pas- 
sage" , problems may be expected to translate the classical results to quantum 
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mechanics. Allcock examined the question some time ago and arrived at a 



negative conclusion 
to be too pessimistic 



but more recently his arguments have been found 
yy many researchers. In particular, his claim on the 
unavoidability of reflection by the detector has been shown to be unfounded 
The very existence of difficulties makes the time of arrival a fascinat- 
ing subject where basic aspects of quantum mechanics have to be considered, 
such as the interpretation of the quantum formalism, the particle-wave du- 
alism, ambiguities in the quantization of a classical variable, the relation 
between classical and quantum mechanics, or the "measurement problem". 
These are all difficult and not completely understood matters so, not sur- 
prisingly, the arrival time, as well as any other problem concerning time 
in quantum mechanics is challenging and, frequently, controversial. Apart 
from the purely conceptual interest there are experiments that provide ar- 
rival times and a theoretical interpretation is required. Indeed, general first 
passage time problems (where the examined variable is not necessarily a po- 
sition) are quite common. An example is the distribution of times required 
for a diatomic molecule to dissociate 0. Several complementary research 
avenues will be discussed in the following sections. They do not necessarily 
conflict with each other, since different results are generally the consequence 
of having posed different questions [O. Work on the various routes is at dif- 
ferent levels of maturity and some ofthe many open questions are indicated. 
Because of length limitations this is not a truly exhaustive review, so we can 
barely touch the surface of many interesting contributions. We hope anyway 
that the uninitiated (and perhaps some experts) will beneflt from the global 
perspective presented here. 



2 Particle trajectories and wave features 

Some interpretations or formulations of quantum mechanics are based on the 
trajectory concept. In the "causal" theory of Bohm [0, or in the stochastic 
interpretation |]14[, the trajectories are supposed to be actual ones and com- 
bine according to standard probability rules. (Instead, in Feynman's path 
integral formalism the contribution of each path is added as a complex am- 
plitude. Other formal paths have also been explored for certain applications, 
for example "Wigner trajectories" , based on a flt of the dynamical equation 
for the Wigner function into the classical Liouville equation, or "Weyl tra- 
jectories" , that retain several useful properties in common with the classical 
ones [l^.) Using one of these theories deflnite, and generally different an- 
swers for characteristic times, such as the traversal time, the dwell time and 
the arrival time are found. The basic problem with these investigations is 
the uncertainty on the ultimate reality of the involved trajectories. Are, for 
example, Bohm trajectories anything more than a possible scenario? There is 
no experimental evidence to support or reject them as actual particle paths. 
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Finding a way to answer this question would be a major breakthrough in the 
foundations of quantum physics. In the "causal" approach to the arrival time 
I J{X, t) I (suitably normalized) is the general arrival time distribution, 
le result holds for arbitrary X even in the presence of a potential barrier. 
Because all trajectory crossings, repeated or not, are taken into account, 
I J(X, t)| does not correspond in general to a first passage time distribution. 

Instead of relying on the problematic "particle" aspect of the quantum 
state, it is also possible to pay attention to times characterizing the evolution 
of the "wave", frequently by means of asymptotic methods. An example 
is the asymptotic "phase time" derived from a stationary phase argument. 
More detailed information is extracted from contour deformations in the 
complex plane and steepest descent methods [|l^. Within this perspective, it 



is natural to define "times of arrival" of special wave features: the main peak, 
the forerunners, the centroid, or a given percentage of probability density. 
These may be very relevant when the detectors are sensitive to them, but 
should not be over-interpreted in terms of "particle paths" . 



3 "Operational" procedures 

The quantum arrival time may also be handled by means of experiments that 
would provide the traversal time classically. Quantum theory should be able 
to predict the statistics of the experimental observations using appropriate 
models for the measurement. (Inversely, it is also possible to derive from the 
statistics of the recorded data the associated quantum-mechanical observable 
[P^.) These "operational methods" may provide bizarre results, possibly 
opposite to classical expectations. In fact different, classically equivalent 
experiments give in general different outcomes. This quantum multiplicity 
associated with a single classical quantity was already emphasized by Bohr 
[ |19|] : Evidence obtained under different experimental conditions cannot he 
comprehended within a single picture, hut must he regarded as complementary, 
in the sense that only the totality of the phenomena exhausts the possible 
information about the objects.''^ 

We shall describe first two operational approaches based on two imple- 
mentations of an "absorbing boundary" . 



3.1 Elimination of norm at a discretized sequence of 
times 

Assume that a wave packet is chopped sharply at a regular sequence of times 
{tj}, i = 1, 2, ... separated b^At in such a way that the part on the right of 



the point X is eliminated 



For a classical ensemble of noninteracting 



particles such a procedure would lead in the limit At — > to the first arrival 
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time distribution since the chopping is effectively acting as an absorbing 
boundary. However, in quantum mechanics the story is quite different. By 
increasing the chopping rate a higher fraction of particles is reflected, and in 
the limit At this process leads to total reflection. 



3.2 Arrival detectors and complex potentials 

The use of complex potentials has been proposed as a way to mimic the 
classical "absorbing boundary" procedure to evaluate first passage times, 
and to model experimental conditions in destructive time of flight experi- 
ments [|TT| , p2| . This approach does not lead to the disappointing result of 



the previous operational method. The processes taking place in the destruc- 
tive detector, such as ionizations, change the structure or internal state of 
the incident particle, i.e., in scattering theory language, they change "the 
channel" . Of course the full measurement can be quite complex and involve 
many degrees of freedom, but a "reduced" Schrodinger equation can always 
be written for the incident channel in terms of a complex (non-hermitian) 
potential, whose precise form is usually modelled phenomenologically. The 
norm of the incident channel, N, is not conserved and may be used to de- 
fine an effective arrival time distribution, proportional to —dN/dt, that will 
in general depend on the particular complex potential (equivalently on the 
detector). Under conditions usually met in time of flight experiments using 
atomic or molecular beams, - detector at asymptotic distance from the scat- 
tering region and particle source - and for good enough absorbers —dN/dt 
can be approximated by J, the quantum mechanical current density with- 
out the absorber. The difference between the first moments of —dN/dt and 
J (the time averages) is the dwell time of the original particle in the com- 
plex potential region |]TT|. In fact the "average" / Jtdt/ J Jdt is an ideal. 



apparatus-independent quantity that can also be obtained in certain condi- 
tions as an average of a time operator (see the discussion below). Regarding 
the flux as the ideal arrival time distribution, for free motion, or at asymp- 
totic distances from scatterers or sources, is an appealing idea because of the 
agreement with the classical expression. However, a quantum mechanical 
state composed by positive momenta is compatible with a negative value of 
J at certain times and positions P,^. Let us first stress that this effect is 
quantitatively negligible for normal practice so that J should be a sat- 



isfactory quantity for the analysis of most experimental data. Nevertheless, 
the fundamental objection to regard J as a true arrival probability remains 
a valid one, and the deviations of experimental distributions from J in the 
quantum backflow regime can and should be studied with the aid of modern 
advances in experimental atomic and optical physics. 

It would be interesting to determine if ideal absorbing conditions (reflec- 
tion coefficients of the complex potential equal to zero for the wave packet 
momentum range) imply a unique potential, or any particular relation be- 
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tween —dN/dt and J (so that an ideal quantity could be defined from the 
operational procedure.) Inverse scattering for real potentials in one dimen- 
sion is a well developed field but very little is known about the "inverse 
scattering problem" for complex potentials. A second open question is the 
explicit construction of models, with additional detector degrees of freedom, 
that justify the phenomenological results on a more fundamental level. 



3.3 Other measurement models 



Recently, Aharonov et al. ||2^ have analyzed a number of simple, ideal- 
ized "toy models" of arrival time measurements. Unlike the previous ap- 
proach, they explicitly include extra "detector" degrees of freedom in the 
(real) Hamiltonian. They conclude that the free particle arrival time cannot 
be measured more accurately than h/ where E is the kinetic energy of the 
particle, but it is necessary to investigate further if this is a model dependent 
result or a fundamental limitation (they argue in favor of the later). 

Schulman has proposed a theory of quantum measurement where the state 
of the studied microscopic system evolves by unitary evolution -including gen- 
erally environment and apparatus in the Hamiltonian- to one of the possible 
eigenstates of the measured observable (corresponding to the result of the 
measurement) |2^. In particular, the particle detection would require the 



localization of the entire particle wave function in the detector [|26[. Thus 
this theory seemingly leads to different results from the ones discussed in 
the previous subsection, where such localization is not assumed. A detailed 
model including apparatus and environment would clarify the actual differ- 
ences further, and the occurrence or not of the localization proposed; also the 
quantitative implications in the calculation of arrival times. Schulman has 
discussed experimental tests that would determine the validity of his theory 



4 Path decomposition expansion and renewal 
equations 

A renewal equation with the form (^) or a related equation in terms of cu- 
mulative probabilities could be formally written in the quantum case and 
solved by Laplace transform by giving some precise meaning to the symbol 
P(x2,a:i,t2 ~ ^i) This meaning however is not at all obvious or clearly 
defined unless some interpretation in terms of trajectories is used where the 
position becomes a Markov process. Localizing the particle around xi, e.g. 
with a Gaussian wave function, is possible P , but the evolution will depend 
on the momentum average and dispersion of the chosen Gaussian. 
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A close quantum relative of the classical renewal equation is the "path 
decomposition expansion" (PDE) of the propagator into a sum over Feynman 
paths classified according to their first passage time at X. For x,x" < X it 
takes the form 

K{x",t"\x,0) = Ko{x",t"\x,0) 

+ f" K(x",t"\X,t')^-^Ko{X,t'\x,0)dt' , (5) 
JO 2m aX 

where Kq is a restricted propagator corresponding to the half-space (—00, X). 
This expression was first derived using Feynman path integrals , and later 
by a more general operator procedure |^ . An even simpler derivation follows 
from the general relation between the propagators K and Kq, corresponding 
to Hamiltonians H and Hq, H = Hq + V , 

K{t) = Ko{t) ~l t K{t - t')VK\t')dt\ (6) 
Ti Jo 

by putting Hq = Q{—x + X)H<d{—x + X) (6 is the Heaviside function) [p9|] . 
By analogy with the renewal equation, it is tempting to consider 

A = J ^^Ko{X,t'\x,0)ij{x,0)dx (7) 

as a "first passage time amplitude" [^,^]. Some caution should however 
be exercised since Kq does not correspond to an absorbing boundary but to 
a reflecting one, so the analogy with (^) is only a partial one. Moreover, 
the squared modulus does not have the correct dimensions nor will generally 
satisfy "probability sum rules" because of interferences between paths taking 
different times. The interferences may however dissapear when coupling the 
particle with an environment. This decoherent effect and the probabilities 
so obtained have been examined by Halliwell and Zafiris |31] within the 
"decoherent histories approach to quantum mechanics" . These authors point 
out that when decoherence is achieved the resultant probabilities depend on 
the mechanism producing decoherence, and insist, quoting Landauer [0,^ , 
that time in quantum mechanics only makes sense if the mechanism by which 
it is measured is fully specified. 

It is illustrative to compare the (appropriately normalized) squared mod- 
ulus of A with the flux J ^9\. For free motion, K and the restricted propa- 



gator Kq are known, and analytical results are available by taking as initial 
state, at time t = 0, a minimum uncertainty Gaussian wave function with 
central position and momentum xq and po, and spatial variance 6"^. If X = 0, 
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where a. = txoh'^. A detailed calculation shows that, putting a = 0, the right 
hand side is proportional, up to time independent factors, to The shapes 
of J and a (normalized) |y4p are close to each other when >> \xQ\h. 



5 Quantization rules and time operators 

The search for "time operators" and the study of their properties has been 
the traditional and most popular approach P,P, p^ - |47[] , even though Pauli 
pointed out the impossibility of a self-adjoint time operator conjugate to 
a Hamiltonian with bounded spectrum |E8|. There are however different 
ways to circumvent this objection by definmg operators which retain at least 
partially the desirable properties of a time observable. Unfortunately, the 
definition of quantum operators associated with classical quantities is not 
justified at present by any fundamental quantization theory, and all known 
quantization rules are essentially heuristic recipes that may provide ambigu- 
ous, non unique, or useless operators for some classical quantities ||4^. It 
is essential in each case to examine the properties of the operators obtained 
and determine their physical content (conditions that they satisfy, domain 
of applicability, and relation to operational procedures and other quantities 
of interest). The connection of the operators with actual measurements is 
frequently obscure: As stated by Wigner |^0|, "There is no rule that would 



tell us which self-adjoint operators are truly observables, nor is there any 
prescription known how the measurements are to be carried out, what appa- 
ratus to use, etc. In a theory with a positivistic undertone, this is a serious 
gap." In principle, any operator resulting from a quantization rule can be 
associated with a property of the state of the system, which may or may 
not be easily measurable or useful. In general, a number of conditions, not 
only motivated by experiments, are imposed to select among the possible 
operators. Claims of uniqueness should then be taken cautiously, since they 
generally reflect the proponent bias towards a group of conditions, which may 
not be satisfactory for certain purposes (and certainly not for the sensibility 
of competitors!) 

Kijowski I^Bj, for the free motion case, imposing a series of conditions 
compatible with the classical arrival time, and limiting the domain to states 
with positive momentum, derived (uniquely within the stated conditions) the 
distribution of arrival times (X = 0) 



n(t;^(o)) 



(m/i)V2 



(9) 



where {p\ip{0)) is the state in momentum representation at t = 0. Note the 
correct behaviour under time translation of the wave function: n[t;?/'(0)] = 



9 



njt — t'; ip{t')]. Related analysis were carried out by Werner and Ludwig 
\ p^ , and more recently by several authors [p|,p9|- i3[1 . This distribution has 
also the remarkable property of having the same nrst moment than the flux 



U{t)tdt = J J{t)tdt , (10) 
and may be formally written as the square | p, where \t) is the eigenstate, 



{P\t) 



mh 



of the time operator p4| , ^p 



t 




(11) 



(12) 



(The symmetrical form —mp ^/"^xp leads in the positive momentum sub- 
space to the same eigenvalue equation and eigenstates as t P,|39[|.) The 
normalization in (pTl) is chosen so that 



{p\t){t\p')dt = 5{p~p'), p,p'>0. 



(13) 



Dealing with states composed by coherent combinations of momenta with ar- 
bitrary sign is more difficult because of the singularity at momentum zero, but 
there should be a theoretical distribution in agreement with an experimental 
arrival time distribution even if the average does not exist. A regularization 
procedure has been proposed I^Of. The results can be very different in this 
case from the ones derived from the Bohm approach and consideration 
of actual experiments in this regime would be of much interest. 



5.1 Phase space quantization techniques 

Quantum states and "observables" can be expressed equivalently in operator 

form (p and G respectively) or by means of various phase space represen- 
tations or images [F{q,p) and g{q,p) respectively] in such a way that the 
expectation value of the operator can be written as 

{G{q,p)) = jj F{q,p)g{q,p)dqdp. (14) 

Sets of four transformations 

p^F, G^g, (15) 
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characterized by a kernel function / can be constructed |]52| - |53| . Each / 
defines the quantization rule g ^ G, 



1 



47r2 



g{q,p)e''^'^''+^P^f{e, T)e'^^^+^P> dqdpdedr , (16) 



and the phase space images of states and observables. The Weyl-Wigner 
formalism, where the state is represented by the Wigner function, , and 
the quantization rule is given by Weyl's prescription, corresponds to / = 1 

g. 

We shall first look at the free motion case. The classical time of arrival 
at point X when the trajectory starts at qQ,p at time t = is given by 



t 



{X - qo)m 



V 



(17) 



Inserting (|T^ in 
Xm 



t 



27r 
m 
2Txi 



the corresponding operators are obtained as 
(P-P)/(0,r)c/pc/r 

e-"^P-^dpdT. 



1 

-e 
p 



fLoiO,r) + f{0,T)[-znT + tq 



1 



(18) 



A wide family of quantization rules, and in particular the ones by Weyl, 
Rivier and Born- Jordan p2| , ^ lead to the same operator we have already 
discussed, expressed now for arbitrary X, 



Xm 
p 



m q 



p 



h 1 
2iP 



Xm 
p 




(19) 



while the standard and antistandard quantizations do not produce a her- 
mitian operator so they will not be discussed further. Several properties 
of this operator are easily proved in phase space using the Weyl-Wigner 
formalism, / = 1. In particular, the relation (p!0| ) can be derived using 
F^{q,p,t) = F^{qQ,p,0) (vahd for free motion, q = qo + tp/m), and not- 
ing that the phase space representative of the fiux operator is the classical 
expression 6{X)p/m. 



6 Other phase space techniques 

Generalizing the free motion arrival time operator is not simple, and only a 
few interaction potentials P6|J77[] or asymptotic distances from a scattering 
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potential have been worked out. The classical expression for the tra- 
jectory may become very involved, and it is rarely explicit (as a function 
of go, Po, and t), so that considering different operator orderings becomes 
cumbersome or impossible in practice. There is the additional difficulty that 
the equation X = q{qo,Po,t) has in general real solution only for a limited 
domain of the initial phase space. 

We shall sketch here how the phase space formalism combined with the 
Heisenberg picture can be used to provide such a generalization, and to deal 
with quantities different from position. The phase space images of the 
Heisenberg operators G(t) depend on the initial phase space point {qo,po) 
and on time, 

p(0)G(t)l = // F{qo,Po,0)g''{qo,Po,t)dqodpo. (20) 



tr 

(In particular the images of the Heisenberg operator for position, x^{qo,po, t), 
are in the classical limit classical trajectories.) Suppose that the equation 
for an arbitrary , 

9''{qo,Po,t) = g, (21) 
where ^ is a predetermined value, has at least one solution for t > 0. We 
can identify the sequence of times t^f\qoiPoiQ)i i =,1,2,... where g^{qo,Poi t) 
"crosses" Q. An average time (note the dependence on /) is then defined for 
the i-th crossing as 

iU) = j^JJ^F{qo,po,0)tf{qo,po,g)dqoPo, (22) 
= JJ^ F{qo,po,0)dqodpo, (23) 

where the domain of integration is restricted to the phase space region Di 
where (^Tj) has an i-th solution, A^j is a normalization constant, and / may 
be tailored in order to satisfy consistency requirements |^ or experimentally 



obtained values. Higher moments can be obtained similarly. 



7 Average "presence" times 

In several of the previous sections the probability flux J has been emphasized 
as an important quantity in relation to the arrival time. We have seen in 
particular that classically it provides the first passage distribution for free 
motion, or for a general case when absorbing boundaries are imposed. The 
probability density p does not play this role. However, average times can 
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also be defined in terms of it. It is appropriate to use a different name for 
tliem, for example average "presence" times, 



/ p{X, t)tdt 
Ip{X,t)dt 



(24) 



They are of interest for detectors sensitive to the presence of the particle 
rather than to the flux. These times, and associated operators have been 



studied by several authors [55-58 



8 Concluding remarks 

Understanding the various aspects of the time of arrival in quantum mechan- 
ics remains an exciting technical and conceptual challenge. Several theoreti- 
cal approaches have been proposed or reviewed. Experiments in non-classical 
regimes (with backfiow, or for motion governed by an interaction potential) 
would provide a much needed reference to refine operational models and 
ascertain the practical relevance of intrinsic quantities or operators. 
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